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ABSTRACT 

Emission-line abundances have been uncertain for more than a decade due to unexplained discrep- 
ancies in the relative intensities of the forbidden lines and weak permitted recombination lines in 
planetary nebulae (PNe) and H II regions. The observed intensities of forbidden and recombination 
lines originating from the same parent ion differ from their theoretical values by factors of more than 
an order of magnitude in some of these nebulae. In this study we observe UV resonance line ab- 
sorption in the central stars of PNe produced by the nebular gas, and from the same ions that emit 
optical forbidden lines. We then compare the derived absorption column densities with the emission 
measures determined from ground-based observations of the nebular forbidden lines. We find for our 
sample of PNe that the collisionally excited forbidden lines yield column densities that are in basic 
agreement with the column densities derived for the same ions from the UV absorption lines. A similar 
comparison involving recombination line column densities produces poorer agreement, although near 
the limits of the formal uncertainties of the analyses. An additional sample of objects with larger 
abundance discrepancy factors will need to be studied before a stronger statement can be made that 
recombination line abundances are not correct. 

Subject headings: planetary nebulae: general — planetary nebulae: individual (He2-138, NGC 246, 
NGC 6543, Tc 1) — ISM: abundances — ultraviolet: ISM 



1. INTRODUCTION 

The analysis of emission line intensities has been used 
to determine nebular abundances for a wide range of 
objects. Standard procedures have been developed in 
which the collisionally excited forbidden lines and high- 
level permitted recombinat ion lines of ions are used 
to determine abundances (jDopita fc Sutherland! 120031 : 
lOsterbrock fc Ferlandll2006h . For most objects heavy el- 
ement abundances are derived from the forbidden lines 
because of their greater strengths compared to the 
fainter recombination lines, which are frequently only 
marginally stronger than the continuum intensity. For 
some of the higher surface brightness gaseous nebulae 
both types of lines have been used to determine the heavy 
element CNO abundances, and they have produced dis- 
crepant results of more than an order of magnitude in 
some objects. 

Each of the two types of lines have certain advantages 
for abundance determinations. The forbidden lines are 
strong, so they are detected from many more ions than 
the recombination lines. The collisional excitation of low- 
lying levels dominates other competing population pro- 
cesses such as fluorescence excitation, charge exchange, 
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and dielectronic recombination. Furthermore, collision 
strengths coupling most of the lower bound levels of ions 
are known to better than 30% accuracy. The largest 
uncertainty in using forbidden line intensities for abun- 
dances is their sensitivity to kinetic temperature that 
results from excitation by electron impact. 

Direct electron recapture populates the higher levels of 
ions and this process has relatively small cross sections. 
Thus, recombination lines tend not to be strong except 
for H and He by virtue of their dominant abundances, 
but they are observable in nebulae from ions of CNO 
and Ne. They have the advantage that recombination 
line intensity ratios are insensitive to temperature and 
density, and the relevant cross sections are believed to 
be known reasonably well. Because recombination cross 
sections are small, however, other excitation processes 
compete with electron recapture in populating the higher 
levels from which these lines are observed. Thus, there 
can be greater uncertainty in the excitation processes 
that are responsible for specific high level permitted lines. 

Electron temperatures and densities are determined 
directly from the relative intensities of forbidden lines 
originating on different levels of the same ion with the 
result that emission spectra have been a major source of 
our knowledge of element abundances of every type of 
emission-line object. The relatively high surface bright- 
nesses of planetary nebulae and a few of the brighter H 
II regions enable the recombination lines of CNO to be 
observed, and in the past decade ion abundances have 
been determined for a number of PNe using both the 
forbidden lines (FL) and the recombination lines (RL) 
from the same ions. Surprisingly, the two types of lines 
have not yielded the same abundances. The differences 
between the RL and FL abundances vary from object to 
object and span the range from 15 percent, i.e., relatively 
good agreement, to factors of 50 and more (|Tsamis et al.l 
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TABLE 1 

Ions with UV Resonance and Optical Forbidden Lines 



Ion 


UV Resonance Transition A(A) a 


Optical Forbidden Line A(A) 


CI 


1277.5*, 1329.6*, 1561.4*, 1657.0* 


4622, 8727, 9850 


p II 


1152.8, 1154.0*, 1301.9, 1310.7*, 1542.3*, 1532.5 


4669, 7876 


SI 


1270.8, 1277.2, 1295.7, 1316.5, 1425.0, 1433.3*, 1474.0, 1483.0*, 


4589, 7725 




1807.3, 1820.3* 




Fell 


1260.5, 1608.5, 1621.7* 


4244, 4359, 5159, 8617 


Ni II 


1317.2, 1370.1, 1454.8, 1709.6, 1741.5 


6667, 7378 


NI 


1199.5, 1200.2, 1200.7 


3467, 5198, 5200, 10398 


OI 


1302.2, 1304.9*, 1306.0** 


5577, 6300, 6364 


SII 


1250.6, 1253.8, 1259.5 


4069, 4076, 6716, 6731, 10320 


S III 


1190.2, 1194.1*, 1201.7* 


3722, 6312, 9069, 9531 



a Single and double asterisks indicate transitions arising from the first and second fine-structure excited level of ground-state 
term, respectively. 



TABLE 2 

Journal of Observations for HST UV Spectroscopy 



Object 


He2-138 


NGC 246 


NGC 6543 


Tc 1 


Central star (V) 


10.9 


11.9 


11.1 


11.4 


Shell surface brightness, S(H/3) 










(erg cm -2 s — 1 arcsec 2 ) 


5.1xl0- 13 


6.2xl0" 16 


8.2xl0" 13 


3.0xl0~ 13 


Diameter of Central Nebular 










Emission (arcsec) 


7 


245 


20 


10 


Radial velocity, heliocentric (km s _1 ) 


-47 


-46 


-66 


-83 


Exposure times (sec): 










1150-1330A 


7x2114 


1967 


2x2150 


2x2072 
8x1365 


1316-1518A 




2x1368 


2440 




1495-1688A 


7x2116 




3x2460 


2077 
4x1367 



2004 iRobertson-Tessi fc Garnettl 12003: iLiu et all 120061: 
Garcia-Roias fc Estebanl 120071) . with the recombination 
line intensities always being stronger than predicted rel- 
ative to the forbidden lines and therefore indicative of 
higher abundances. These discrepancies have been the 
subject of many studies which have given rise to a large 
literature on the subject, but they are still not un- 
derstood. Until they are resolved some doubt is cast 
upon the normal methods by which collisionally excited 
forbidden lines are used to derive element abundances. 
The differences cannot be due to incorrect atomic data 
since this would cause the magnitude of the discrepan- 
cies to be roughly the same for all objects. Current res- 
olutions to the discrepancies problem have focused on 
temperature fluctuat ions in the nebulae (|PeimbertJll967l : 
iPeimbert et al . 2004) a nd dense inclusions that are hy- 
drogen deficient (|Liu et al"1l2000D . 

A completely independent, alternative method of ob- 
taining abundances for nebular gas does exist and can 
be used as an independent check of emission-line abun- 
dances. It involves observing the absorption lines pro- 
duced by the foreground nebular gas in the spectrum 
of an embedded or background star to determine column 
densities. Most of the absorption lines occur in the ultra- 
violet because low density gas occupies the ground state 
and the resonance lines of the most cosmically abundant 
ions fall in the UV. Thus, a space telescope with a high 
resolution spectrograph is required to study these ab- 
sorption lines with sufficient resolution to yield reliable 
column densities. 



IPwa. Mo. fc Pottaschl <|1984D and iPwa. Pottasch. Mol 
(|1986l ) made the first attempts to obtain ion abundances 
in planetary nebulae by this method, using the IUE high 
resolution spectrograph (R=15,000) to measure absorp- 
tion line equivalent widths from which column densities 
were determined. For the two PNe whose central stars 
were bright enough to be studied with IUE, the relatively 
few unsaturated lines for which they were able to obtain 
column densities belonged to ions that do not have de- 
tectable emission lines in the optical. Hence, although 
they did determine relative abundances from absorption 
line data, they were unable to make a direct comparison 
between independently derived absorption and emission 
line abundances for the same ions. 

The high resolution spectrographs of HST increase the 
number of nebulae for which reliable column densities 
can be obtained, and offer the possibility of resolving 
the question of the abundance discrepancies. The more 
abundant heavy element ions having resonance lines in 
the UV between AA1150-1800A accessible with the STIS 
spectrograph and also having detectable forbidden lines 
at optical wavelengths are listed in Table [TJ The abun- 
dances of these ions can be determined independently 
from ground and space telescopes by completely differ- 
ent methods and then compared, albeit having separate 
sight lines and path lengths, e.g., the line of sight to the 
star does not probe the rear part of the nebula. 

Specifically, the column densities of individual ions can 
be found from the UV absorption spectra, while emission 
measures are derived from the forbidden and permitted 
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Fig. 1. — A montage of UV absorption line profiles from 
HST/STIS spectra of the central star of He2-138. The foreground 
ISM absorption is centered around velocity = -10 km s — 1 , and the 
PN shell absorption is centered around velocity = -60 km s — 1 , as 
indicated by the vertical dashed line. 

nebular emission lines. For each ion the column density 
and emission measure differ only by the multiplicative 
factor of the electron density in the emission measure. 
Since standard nebular diagnostics provide a direct de- 
termination of the density appropriate for each ion de- 
pending upon its ionization level, a direct comparison can 
be made between the absorption column density of the 
ion and its emission measure as derived from the different 
emission lines. This procedure should demonstrate which 
emission lines, forbidden or permitted, yield abundances 
most consistent with those from the UV absorption lines. 

An initial study of abundances determined from UV 
absorption lines in the central star vs. those found from 
nebular emissi on line intensit i es for the PN IC 418 was 
attempted by IWilliams et al . (2003). For the four ions 
S +2 , S + , Ni + , and Fc + , and netural oxygen 0°, for which 
relative abundances could be determined independently 
from both methods, rough agreement was found. How- 
ever, the uncertainties were too large for meaningful con- 
clusions to be drawn. 

We report here on an observing program which at- 
tempts to resolve the discrepancies between the forbid- 
den and permitted emission line intensities by making a 
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Fig. 2. — A montage of UV absorption line profiles from the 
star of NGC 246 from the HST/STIS data. The ISM absorption is 
centered near velocity = 0, and the PN shell absorption is centered 
around velocity = -80 km s —1 , as indicated by the vertical dashed 
line. The broad N V absorption trough is indicative of an outflow- 
ing wind from the central star, which has a heliocentric velocity of 
v* = -46 km s _1 . 

UV absorption line analysis that independently serves to 
validate emission line results. We have obtained high 
resolution UV spectra of four PNe central stars with 
HST/STIS, and visible spectra of three of the associ- 
ated nebular shells from Las Campanas and KPNO. The 
UV observations and absorption analysis are described 
in Sj3J and the optical emission spectra and analysis are 
presented in Sj4j The relative column densities from the 
two methods are compared and interpreted in |J5] 

2. OBJECT SAMPLE AND OBSERVING PROGRAM 

Column densities determined from absorption lines 
are most reliable when the lines are well-resolved and 
have ample signal-to-noise to define the continuum, thus 
brighter central stars are advantageous. Absorption 
lines originating in nebular gas are frequently seriously 
blended with and obliterated by stronger absorption from 
the same transitions caused by intervening ISM gas along 
the same line of sight. Unambiguous measurement of ab- 
sorption from the nebular gas therefore requires a neb- 
ular radial velocity differing by at least 50 km s _1 from 
that of the Local Standard of Rest to shift the nebular 
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Fig. 3. — Same as Figure[T]except for the central star NGC 6543. 
The three vertical dashed lines represent distinct absorption com- 
ponents from the nebular shell. 

absorption out of the corresponding stronger ISM com- 
ponent. Optimal candidates for emission study are pref- 
erentially high surface brightness objects, thus favoring 
PNe over the lower surface brightness H II regions. It 
would be advantageous to include in our sample some 
PNe for which the largest FL and RL abundance dis- 
crepancies have been determined, however the few PNe 
that have been established to have differences of more 
than a factor of ten either have (a) central stars that are 
too faint in the UV, (b) very low surface brightnesses, 
or (c) radial velocities that are not sufficiently different 
from the LSR to avoid confusion between the nebular 
shell and ISM absorption lines. 

The sample of known PN e satisfying the optim al cri- 
teria for study is given in IWilliams et ail (|2003f) . and 
is not large. We identified four PNe that satisfy these 
criteria and which seem well suited for a combined UV- 
visible study, viz., He2-138, NGC 246, NGC 6543, & 
Tc 1, whose central stars are sufficiently bright that UV 
observations with HST at high spectral resolution would 
produce acceptable spectra in reasonable exposure times. 
We acquired spectra of the central stars of these PNe in 
the UV with HST/STIS and then subsequently observed 
the nebular shells along adjacent sight lines in the visible 
with ground-based telescopes to obtain line intensities for 
the emission-line analysis. 

3. ABSORPTION-LINE ANALYSIS 

3.1. STIS Observations 

HST/STIS was used to obtain spectra of the central 
stars of He2-138, NGC 246, NGC 6543, and Tc 1 in the 
high resolution mode, i.e., grating E140H with a resolu- 
tion of 3 km s , in three separate settings that covered 
the wavelength region 1150-1690 A. Exposure times that 
produced a continuum signal-to-noise level of S/N«15 
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Fig. 4. — Same as Figure [T] except for the central star of Tc 1. 

over the entire wavelength regime for each grating set- 
ting were adopted. The observations were made in 2005 
(Cycle 12), and the relevant properties of our targets 
and the journal of the STIS observations is given in Ta- 
ble [U Regrettably, STIS failed and became inoperative 
before our observing program could be completed, thus 
we did not succeed in executing all of our planned obser- 
vations. Only partial data exist for each of the central 
stars. The spectra were reduced using the mos t recent 
versio n of CALSTIS procedures and algorithms (|Lindlerl 
1998), and a montage of resonance line profiles from the 
final reduced spectra that include ions for which we also 
subsequently observed nebular forbidden emission lines 
is shown in Figures [TJ3] for the four PNe. 

3.2. UV Line Measurements 

Table [3] lists the absorption lines that we measured in 
the central star spectra of the four planetary nebulae. 
Following the name of each target in the subheaders of 
the table, we list the radial velocity of the central star 
and the velocity interval that covers the strongest ab- 
sorption features that we identify as arising from the 
planetary nebula shell. Weaker features often spanned 
smaller velocity intervals and the measurements of these 
lines, along with those undetected, were taken over the 
more restricted ranges. 

We defined the continuum levels by fitting Legendrc 
polynomials to the fluxes on either side of eac h line, using 
the methods refined by ISembach &: Savage] (|1992D . In 
Figure [5] we present a portion of the UV spectrum of the 
He2-138 central star which shows the final continuum 
fit, together with the envelope defined by la excursions 
from the fit, used in the determination of the absorption 
intensities in the apparent optical depth analysis of the 
Ni II A1317.217 line. 

For each absorption feature, within the errors of the 
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Fig. 5. — A segment of the STIS spectrum that covers the Ni II 
line at 1317.217 A for the central star in He2-138. The large feature 
centered at 1317.15 A arises from foreground interstellar material. 
The small feature centered 1316.96 A is the one that is relevant to 
our study, since it arises from the nebula. The adopted continuum 
level is shown by a curved line with a cross-hatched overlay that 
shows the lcr uncertainty in its placement. 

continuum fitting there is an acceptable range for the re- 
constructed intensity levels, and limits for this range de- 
fined the errors in the line measurements attributable to 
continuum uncertainties. For both the equivalent width 
measurements and the evaluations of column densities 
using the apparent optical depth (AOD) method (i j3.3p . 
we combined these continuum uncertainties in quadra- 
ture with the uncertainties due to photon-counting noise 
to arrive at a net error of the quantity being measured. 
In some instances, lines could be measured twice because 
they appeared in two adjacent echelle orders. 

While our principal objective was to obtain column 
densities for ions in the nebular shells for comparison 
with emission line strengths from the same ions, we nev- 
ertheless have included in Table [3] measurements for ab- 
sorption lines from ions for which there were no emission 
lines detected in our ground-based spectra fi j4.2[) . We feel 
that it is prudent to include these lines for the benefit of 
future more general studies of the relative abundances of 
atoms and ions in the planetary nebulae shells. 

Uncertain column densities result from either very 
weak lines that are marginally detected or strongly sat- 
urated lines. A few of the ions have multiple lines with a 
range of f-values that provide for reliable column densi- 
ties. We avoided lines that were so badly saturated that 
their resulting lower limits for the column densities would 
be so much lower than the actual values that they would 
have no real value for any study. For some species, e.g., 
Mg II, only the strongest line appeared above the noise; 
for such cases we could not measure the weaker line. Us- 
ing the same methods as for lines that were visible in 
our spectra, we evaluated intensity upper limits within 
the wavelength intervals where certain lines of interest 
might be expected, but which were either marginally de- 
tectable or not visible at all. Sometimes these measure- 
ments yielded negative equivalent widths, although with 
magnitudes comparable to or much less than the errors, 
and these determinations are ultimately useful in provid- 
ing upper limits for the column densities (see footnote b 
to Table El). 

3.3. Absorption Column Densities 



Column densities N were derived by integrating 
over velocity the apparent optical depths r a (v) — 
^■[Icontiv) / I{v)\ and evaluating the quantity 

N = J mdt = m e c/(ne 2 f\) J T a {v)dv , (1) 

where n± is the density of the ion in the lower level, and 
the numerical value for the expressio n in front of the inte- 
gral is 3.77 x 10 14 c m" 2 (km s -1 ) -1 (Sav age fc Sembachl 
ll99lHJenkini ll996). The results for all of the ions with 
reliable determinations are given in Tabled and it should 
be emphasized that the column densities refer only to 
those ions that occupy the lower level of the transition. 
Errors in the column densities may arise from three dif- 
ferent sources: (a) photon-counting noise, (b) errors in 
the definition of the continuum level, and (c) errors in 
the adopted zero intensity level. If the random devia- 
tions of intensity arising from statistical fluctuations in 
photon counts are expressed as o~i( v ) , an approximation 
for the error in T a (v) is simply 

a r a (v) = <7l(v)Icont(v)/I(v) , (2) 

which is reasonably accurate as long as the quantity is 
much less than unity. iJenkins fc Trippl (|2001h found that 
for a signal-to-noise ratio of about 20 at the continuum 
(which applies to nearly all of our spectral lines) and a 
Gaussian error distribution, the approximation expressed 
in equation[2]is good as long as I{v)/I cont {v) > 0.15. We 
have indicated which lines appearing in Table [3] violate 
this condition at the maximum level of absorption. For 
these cases, the upper error bounds may need to be in- 
creased to somewhat larger values than those listed. 

The errors in optical depth cr Ta r v \ that arise from pho- 
ton counting are uncorrelated from one spectral element 
to the next, while the systematic error arising from a 
misplacement of the continuum is an effect that is usu- 
ally coherent over the extent of an absorption feature. 
For this reason, the noise errors for successive spectral 
elements were added together in quadrature before they 
were combined as a group with the global uncertainty in 
line strength caused by inaccuracies in the definition of 
the continuum level. Since the errors arising from photon 
counting and continuum misplacement are uncorrelated, 
it is appropriate to add them together in quadrature. At 
the bottoms of severely saturated lines, we found that 
the intensities deviate from zero by less than 1% of the 
continuum intensity. Thus, anomalies arising from errors 
in zero corrections are insignificant compared to short- 
comings of the approximation in equation O 

Except for some strong lines of Si II and Si II* in the 
spectrum of NGC 6543, the values of TV obtained for lines 
of different strength generally agreed with each other. In 
this exception, the fact that the stronger lines yielded 
lower column densities than the weaker ones for these 
two ground state levels indicates that there are some un- 
resolved saturated absorptions that make all of the eval- 
uations of N using equat ion Q] underest imate slightly the 
true value of N for Si II penkinslll99ll . 

4. EMISSION-LINE ANALYSIS 

4.1. Optical Spectroscopy 

The four PNe studied here were observed in the visi- 
ble with ground-based telescopes to measure intensities 
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TABLE 3 

Equivalent Widths and Column Densities 

A a log/A a Species W x b logJV b 

(A) (mA) (cm- 2 ) 





He2-138 


(u» = —47 km s 1 ; 


- 73 < tip N < -40 km 


s- 1 ) 






9 3Q9 

z.oyz 


r 1 t 


98 -t- A A 
Zo.u ZH <±.*± 


1 9 Q7 


^tu.ui , — yj.vjo) 




2.082 


C I 


2.6 ± 1.8 


12.25 


(4-0 22 —0 47) 


1277.245 


2.037 


C I 


11.2 ± 1.8 


13.00 


(+0.07, -0.08) 


1548.204 


2.468 


C IV 


32.5 ± 4.2 


12.99 


(+0.06, -0.07) 


1550.781 


2.167 


C IV 


15.7 ± 4.0 


12.93 


(+0.10, —0.13) 


1550.781 


2.167 


C IV 


11.8 ± 2.6 


12.79 


(+0.09, -0.11) 


1304.858 


1.795 


O I* 


75.1 ± 1.2 


14.49 


+0.05, -0.05) c 


1306.029 


1.795 


O I** 


64.9 ± 1.4 


14.38 


+0.04, -0.04) c 


1306.029 


1.795 


O I** 


67.7 ±2.2 


14.20 


+0.08, -0.10) c 


1239.925 


-0.106 


Mg II 


5.0 ± 1.3 


14.79 


(+0.10, -0.14) 


1239.925 


-0.106 


Mg II 


4.2 ± 1.1 


14.72 


(+0.10, -0.14) 


1670.789 


3.463 


Al II 


91.4 ± 20.5 


12.69 ( 


f 0.12, — 0.17) c,d 


lOUI.O ( \J 


2.052 


Si II 


QQ Q _)_ 1 n 


14.57 ( 


■uc\ ozl — n n^ c ' d 


1309.276° 


2.052 


Si II* 


69.6 ±2.3 


13.98 


[+0.02, -0.02) f 


1309.276 6 


2.052 


Si II* 


71.5 ±2.4 


13.98 


;+0.02, -0.02) f 


1152.818 


2.451 


P II 


31.2 ±3.8 


13.28 


;+o.o9, -o.n) f 


1301.874 


1.219 


P II 


4.1 ± 1.1 


13.36 


(+0.11, -0.14) 


1153.995 


2.331 


p JJ** 


11.2 ±4.5 


12.83 


(+0.14, -0.22) 


1295.653 


2.052 


S I 


1.6 ±1.6 




< 12.58 


1250.578 


0.832 


S II 


74.9 ± 2.0 


15.49 


+0.08, -0.09) c 
+0.10, -0.14) c 


1250.578 


0.832 


S II 


71.7 ± 1.2 


15.42 


1190.203 


1.449 


S III 


73.6 ± 2.3 


14.76 


;+0.12, -0.16) f 


1190.203 


1.449 


S III 


70.0 ± 1.6 


14.70 


;+0.02, -0.02) f 


1201.729 


0.626 


S III** 


16.3 ±3.2 


14.63 


(+0.08, -0.10) 


1197.184 


2.414 


Mn II 


9.4 ±2.5 


12.66 


(+0.11, -0.15) 


1197.184 


2.414 


Mn II 


10.8 ±2.2 


12.68 


(+0.08, -0.10) 


1199.391 


2.308 


Mn II 


6.2 ± 2.5 


12.53 


(+0.14, -0.22) 


1608.451 


1.968 


Fe II 


55.0 ±2.3 


13.95 


;+0.03, -0.04) f 


1317.217 


1.876 


Ni II 


4.1 ± 1.3 


12.70 


(+0.12, -0.16) 


1317.217 


1.876 


Ni II 


5.3 ± 1.1 


12.81 


(+0.09, -0.11) 


1237.059 


3.183 


Ge II 


9.7 ± 1.3 


11.86 


(+0.06, -0.06) 


1237.059 


3.183 


Ge II 


11.6 ± 1.2 


11.94 


(+0.05, -0.05) 


1235.838 


2.402 


Kr I 


3.4 ± 1.9 


12.14 


(+0.18, -0.33) 


1235.838 


2.402 


Kr I 


4.6 ± 1.2 


12.25 


(+0.10, -0.14) 





NGC 246 


(u„ = —46 km s 1 ; 


- 95 < ti PN < -60 km f 




1277.245 


2.037 


C I 


0.2 ± 1.0 


< 12.29 


1334.532 


2.234 


C II 


12.9 ±1.1 


12.84 (+0.04, -0.04) 


1334.532 


2.234 


C II 


14.3 ± 1.0 


12.89 (+0.03, -0.03) 


1335.708 


2.234 


C II* 


15.8 ± 1.0 


12.93 (+0.03, -0.03) 


1199.550 


2.199 


N I 


-0.7 ± 1.0 


< 12.04 


1302.169 


1.796 


O I 


4.3 ±2.1 


12.80 (+0.17, -0.28) g 


1302.169 


1.796 


O I 


0.6 ± 1.7 


< 12.78 g 


1304.858 


1.795 


O I* 


-1.7 ± 1.3 


< 12.44 


1260.422 


3.171 


Si II 


6.5 ± 2.6 


11.61 (+0.14, -0.22) 


1264.738 


3.125 


Sill* 


-0.3 ± 1.0 


< 11.14 


1206.500 


3.293 


Si III 


49.7 ±3.5 


12.52 (+0.03, -0.03) 


1393.760 


2.854 


Si IV 


145.1 ±2.3 


13.50 (+0.01, -0.01) f 


1402.773 


2.552 


Si IV 


97.6 ± 1.3 


13.51 (+0.01, -0.01) 


1259.518 


1.320 


S II 


3.5 ±2.2 


13.19 (+0.20, -0.40) 


1190.203 


1.449 


S III 


3.4 ± 1.3 


13.08 (+0.14, -0.21) 


1194.058 


1.325 


S III* 


-1.7±2.5 


< 13.32 


1317.217 


1.876 


Ni II 


2.1 ± 1.8 


12.40 (+0.26, -0.71) 


1370.132 


1.906 


Ni II 


-0.6 ± 1.9 


< 12.61 



NGC 6543 (v, = -66km s" 1 ; - 103 < v PN < -74 km s" 1 ) 



1277.245 


2.037 


C I 


1.1 ± 1.4 


< 12.52 


1200.223 


2.018 


N I 


2.0 ± 1.7 


12.29 (+0.25, -0.69) 


1304.858 


1.795 


O I* 


1.7 ±1.3 


12.39 (+0.25, -0.65) 


1306.029 


1.795 


O I** 


-0.3 ± 1.8 


< 12.73 


1306.029 


1.795 


O I** 


-1.4 ± 1.3 


< 12.47 


1670.789 


3.463 


Al II 


6.5 ±4.5 


11.22 (+0.22, -0.48) 
13.11 (+0.02, -0.02) c 


1260.422 


3.171 


Si II 


85.9 ± 1.6 


1526.707 


2.307 


Si II 


35.7 ± 2.0 


13.23 (+0.03, -0.03) 


1304.370 


2.052 


Si II 


19.8 ± 1.4 


13.25 (+0.03, -0.03) 
13.17 (+0.02, -0.02) c 


1264.738 


3.125 


Si II* 


88.8 ±0.9 
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TABLE 3 

(continued) 





log/A a 


Species 


W x b 


log JV b 


(A) 






(mA) 


(cm- 2 ) 






NGC 6543 


(d* = —66 km s ; 


- 103 < D PN < 


-74 km s" 1 ) 


1533.432 


2.307 


Si II* 


43.5 ±2.9 


13.33 (+0.03, -0.03) 


1265.002 


2.171 


Si II* 


33.2 ± 1.5 


13.39 (+0.02, -0.02) 


1309.276 


2.052 


Si II* 


22.6 ± 1.0 


13.31 (+0.02, -0.02) 


1152.818 


2.451 


P II 


5.8 ±6.2 


< 12.81 


1259.518 


1.320 


S II 


9.6 ±2.4 


13.65 (+0.10, -0.12) 


1253.805 


1.136 


S II 


6.3 ± 1.3 


13.64 (+0.08, -0.10) 
14.92 (+0.29, -1.33) c 


1190.203 


1.449 


S III 


93.2 ± 1.4 


1194.058 


1.325 


S III* 


87.2 ± 1.7 


14.96 (+0.02, -0.03) f 


1194.058 


1.325 


S III* 


86.5 ± 1.6 


14.95 (+0.02, -0.02) f 


1197.184 


2.414 


Mn II 


1.3 ±2.1 


< 12.32 


1199.391 


2.308 


Mn II 


0.2 ± 2.2 


< 12.38 


1608.451 


1.968 


Fe II 


0.1 ±4.0 


< 12.84 


1317.217 


1.876 


Ni II 


-0.7 ± 1.6 


< 12.56 


1370.132 


1.906 


Ni II 


-1.2 ± 2.2 


< 12.64 



Tel (v* = -83km s" 1 ; - 128 < v PN < -86 km s" 1 ) 



1560.309 


2.082 


C I 


-3.2 + 4.1 


< 12.65 


1277.245 


2.037 


C 1 


-1.9 ±2.4 


< 12.54 


1199.550 


2.199 


N I 


-2.1 ±3.1 


< 12.53 


1302.169 


1.796 


O I 


27.1 ± 2.4 


13.63 (+0.04,-0.04) h 


1306.029 


1.795 


O I** 


0.5 ± 2.0 


< 12.83 


1239.925 


-0.106 


Mg II 


0.2 ±3.2 


< 14.94 


1239.925 


-0.106 


Mg II 


7.2 ±3.8 


14.94 (+0.18, -0.32) 


1670.787 


3.463 


Al II 


74.8 ± 22.9 


12.34 (+0.12, -0.16) 


1260.422 


3.171 


Si II 


110.7 ±2.1 


13.22 (+0.02, -0.02) c ' d 


1193.290 


2.842 


Si II 


66.9 ±3.7 


13.17 (+0.03, -0.03) 


1190.416 


2.541 


Si II 


38.1 ±4.4 


13.15 (+0.05, -0.05) 


1304.370 


2.052 


Si II 


22.3 ± 1.2 


13.29 (+0.02, -0.03) 


1264.738 


3.125 


Si II* 


99.8 ± 1.7 


13.18 (+0.02, -0.02) c 


1265.002 


2.171 


Si II* 


22.9 ±2.4 


13.20 (+0.04, -0.05) 


1309.276 


2.052 


Si II* 


22.8 ± 1.7 


13.29 (+0.03, -0.03) 


1152.818 


2.451 


P II 


7.3 ±6.9 


12.52 (+0.25, -0.66) 


1259.518 


1.320 


S II 


10.3 ±2.0 


13.67 (+0.08, -0.09) 


1190.203 


1.449 


S III 


79.9 ±2.6 


14.65 (+0.02, -0.02) 


1194.058 


1.325 


S III* 


49.5 ±5.6 


14.46 (+0.05, -0.06) 


1194.058 


1.325 


S III* 


48.4 ±6.2 


14.46 (+0.06, -0.07) 


1197.184 


2.414 


Mn II 


25.8 ±4.8 


13.06 (+0.08, -0.09) 


1197.184 


2.414 


Mn II 


25.9 ±2.9 


13.05 (+0.05, -0.05) 


1608.451 


1.968 


Fe II 


-2.9 ±8.5 


< 13.12 


1317.217 


1.876 


Ni II 


0.4 ± 1.4 


< 12.60 


1237.059 


3.183 


Ge II 


-5.3 ±3.2 


< 11.43 


1237.059 


3.183 


Ge II 


-2.4 ±3.8 


< 11.63 


1235.838 


2.402 


Kr I 


9.0 ±2.6 


12.53 (+0.11, -0.15) 



a Wavelengths and line strengths from [Morton (2000, 2003), exc ept fo r the /-values of Ni II, for 
which we have adopted the values measured by Jenkins & Tripp (2001). Transitions for individual 
species are arranged according to decreasing line strength. This was done in order to make it easy 
to identify trends (strong lines indicating smaller N than weak ones) that signify possible unresolved 
saturated components that could lead to underestimates of column density using the AOD method 
(Savage & Scmbach 1991; Jenkins 1996). Duplicate entries signify independent measurements made in 
different echelle orders. 

b Listed errors represent ±ler deviations and include uncertainties caused by both photoevent statistical 
fluctuations and continuum uncertainties, combined in quadrature. When a measurement of W\ yields 
a value that is below the calculated la error in W\, we state the formal measurement of W\ and its 
error, but then we fo l low w ith an evaluation of a 2er upper confidence bound for the real W\ using the 
method of Marshall (1992) for interpreting marginal detections (or nondetections) of quantities that 
are known not to ever be negative. The listed upper limit for N is calculated from this W\ limit using 
the formula for weak lines (i.e., the linear part of the curve of growth). 

c The line is strongly saturated (central optical depth tq > 4), but without a flat bottom that would 
signify severe saturation. The formal errors listed here may not accurately reflect the true errors. In 
our apparent optical depth integrations, occasional deviations in r a (v) that exceeded 5.0 were simply 
set to equal to 5.0. 

d The right-hand portion of the profile is partly blended with the left-hand portion of the absorption 
arising from foreground material in the general interstellar medium. Thus, the errors could be somewhat 
larger than those derived formally (and stated here). 

c Si II* was also recorded at 1533.4 A, but this line is strong enough to have a small portion of its profile 
strongly saturated. Since our recording of the 1309.3 A feature is of excellent quality (and appeared in 
two orders), we decided not to measure the stronger line. 
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TABLE 3 

(continued) 



f At the bottom of the line, the intensity relative to the continuum is less than 0.15. For our representa- 
tive S/N = 20 (at the continuum), the approximation given in equation [2] starts to become inaccurate. 
For this reason, the upper bound for log N should be increased slightly beyond the value listed here. 
g Absorption by the 1301.874 A transition of P II caused by foreground gas interferes with the O I 
feature. However, we could compensate for this by dividing the spectrum by the P II profile at 
1152.818 A after its strength had been reduced to reflect the fact that the 1301.874 A line is weaker. 
(All intensities in the strong profile were taken to a power equal to the ratio of the lines' values of /A.) 
h The right-hand portion of this profile is partly blended with the left-hand side of an absorption arising 
from the 1301.874 A transition of P II created by the foreground gas. 
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Fig. 6. — Images of program PNe showing locations, orientations, and size of slits for (a) He2-138, (b) Tel, and (c) NGC 6543. In all 
figures North is up and East is to left. Panels (a) and (c) are HST images. The shell emission sampled by our spectra would be represented 
by these images convolved with a 1—1.3 arcsec PSF. 



of the nebular emission lines. The lower surface bright- 
nesses of the nebular shells compared with the central 
stars dictated that the nebular spectra sight lines be po- 
sitioned no closer than 2 arcsec from the central star in 
order to avoid unacceptable levels of contamination by 
scattered light from the brighter star. However, our goal 
was to estimate the intensity directly along the line of 
sight to the central star since this is the position of the 
gas that produces the absorption lines. Our approach 
was to take spectra in two positions symmetrically placed 
on either side of the central star and to then average 
together the two spectra after they had been flux cali- 
brated. The flux average serves to compensate for sur- 
face brightness fluctuations, but does not represent the 
flux along the sight line to the central star if there is a ra- 
dial gradient in surface brightness away from the central 
star due to the three-dimensional structure of the shell. 

Initially, a reconnaissance was carried out in 2004 at 
low spectral resolution (300-500 km s _1 FWHM) using 
the Gold Spectrograph on the 2.1m telescope at Kitt 
Peak National Observatory (KPNO) to observe NGC 246 
and NGC 6543, and the Wide Field CCD (WFCCD) 
camera on the du Pont 2.5m telescope at Las Campanas 
Observatory (LCO) to observe He2-138, NGC 246, and 
Tc 1. We found that the emission lines in NGC 246 are 
much too faint for accurate spectrophotometric measure- 
ments of any of the weaker emission lines; only a few of 
the very strongest lines such as [O III] A5007 were vis- 
ible even in long exposures. This object was therefore 
removed from our program, but we provide the UV in- 
formation from the central star spectrum in this paper 
because of its potential use for abundance studies. The 
visible spectra of the nebulae showed that line blend- 
ing dictated the need for much higher resolution data to 



show many of the weaker emission lines of interest in the 
other three objects. This led us to obtain 15-20 km s _1 
FWHM resolution echelle spectra in 2005 with the echelle 
spectrographs on the 4m Mayall Telescope at KPNO and 
the 2.5m du Pont Telescope at LCO. 

4.1.1. Observations at LCO 

The LCO echelle spectrograph uses a prism as a 
cross disperser and covers the wavelength range AA3480- 
10,150A in 64 orders, with increasingly large wavelength 
gaps between orders beyond A8000A. On each of the two 
nights 8, 9 June 2005 UT, we observed He2-138 and Tc 
1 in two different slit positions symmetrically placed on 
either side of the central star, using a 2 arcsec wide x 
4 arcsec long slit and offsetting at an angle so that the 
slit would include the brightest part of the nebulae. The 
slit orientations were at right angles to the directions of 
the offsets. The two slit positions for Tc 1 were 1.0 arc- 
sec S, 2.7 arcsec W of the central star, and 1.0 arcsec N, 
2.7 arcsec E. For He2-138, the two slit positions were 2.7 
arcsec N, 1.25 arcsec E and 2.7 arcsec S, 1.25 arcsec W 
of the central star. The slit positions used for our ob- 
servations are shown in Figure [5J overlayed on the best 
available images that we could find for these objects. 

For Tc 1 the spectra were taken well inside the outer 
edge of the nebula. However, for the smaller He2-138 
we were not able to simultaneously avoid the scattering 
effects of the central star and sample the shell well inside 
its outer edge, so the slit had to be placed near the outer 
edge of the nebula. Figure 6a shows a Hubble Telescope 
image of He2-138 with the spectrograph slit positions su- 
perposed. Atmospheric seeing resulted in shell emission 
filling the central portion of the slit. We detected flux 
along the central 3 arcsec of the slit, but variations in 
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Fig. 7. — Nebular emission line profiles from the LCO echelle 
spectrum of He2-138. The panels show the forbidden emission 
lines listed in Table [7] that were used to determine forbidden line 
abundances. The corresponding lines are shown at zero velocity in 
each panel. The linear vertical scale is in units of erg cm -2 s — 1 
A -1 with the bottom abscissa of each panel corresponding to zero 
flux and the top corresponding to the value of F^ printed inside 
the panel. 
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Fig. 8. — Emission line profiles from the KPNO echelle spec- 
trum of NGC 6543. The panels and scaling are as in Figure [7] 
except that the right-hand column shows four of the stronger O II 
recombination lines from Table ITOl The stronger red [O I] emission 
components at +70 km s" 1 are atmospheric airglow lines. 

tracking and seeing caused the intensity to vary among 
the exposures. These variations make our determinations 
of the emission measure for this object considerably less 
certain than for the other two PNc. The spectral res- 
olution was 15 km s^ 1 FWHM over most of the range, 
but degraded to 20 km s _1 at the extreme ends. We 
extracted spectra from each position, and after applying 
the proper flux calibration averaged the two extracted 
spectra together. We added together seven 1200 sec ex- 
posures at each slit position over two nights to measure 
the weak lines, and used pairs of 30 or 60 sec exposures at 
each position to measure the strong emission lines which 
would otherwise be saturated. 

Only the second night was photometric so we used ob- 
servations of_th£_s_tandarf_stars HR 4468, HR 4963 and 
HR 5501 (jHamuv et alJll994ft made through a 8x8 arc- 
sec slit on the second night to calibrate both nights. We 
applied this calibration in three steps. We first added 
together the raw counts from all of the different long 
exposures, which is the optimal weighting for detecting 
weak lines, and then flux calibrated that spectrum using 
a mean airmass. This gave a high signal-to-noise ratio 
spectrum that was affected by a wavelength-independent 
attenuation and by a slight error in the wavelength de- 
pendence of the extinction correction. The second step 
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Fig. 9. — Emission line profiles from the LCO echelle spectrum 
of Tc 1. The panels and scaling are as in Figure [jj with the narrow 
[O I] components at +100 km s _1 again due to atmospheric [O I] 
emisison. 

was to properly flux calibrate a single long exposure in 
each slit position on the photometric night using the cor- 
rect airmass. Finally, we measured the fluxes of the same 
intermediate strength emission lines in both spectra, and 
used the flux ratios to correct the high S /N spectrum to 
match the flux scale of the single well-calibrated spec- 
trum. This procedure provides the best calibration for 
our observing circumstances and results in absolute spec- 
trophotometry accurate to better than 8 percent over the 
whole wavelength range. 

4.1.2. Observations at KPNO 

The KPNO echelle spectra of NGC 6543 were taken 
over the three nights 18-20 June 2005 UT. We used 
the UV camera on the 4m Mayall Telescope Cassegrain 
echelle spectrograph with echelle grating 79-63 and cross 
disperser 226-1 with two different setups, each giving 20 
km s _1 resolution. The blue setup, used for the first two 
nights, covered the wavelength range AA3200-5300A. It 
used the cross disperser in second order with a CUSO4 
order separating filter. We then switched to a red setup 
covering the range AA4750-9900A in first order of the 
cross-disperser grating, using a GG495 order separating 
filter. For calibration we observed stand ard stars HR 
4468, HR 4963, HR 5501 and HR 8634 from lHamuvetai] 
(1994), measured through a 6 arcsec wide slit. The nights 
with the blue setup were photometric, while there were 
some clouds present when we used the red setup. We 
used the measured strengths of emission lines in the over- 
lapping sections of the red and blue spectra to scale the 
fluxes for the red spectra to match those of the blue spec- 
trum. 

As was done at LCO, spectra were taken in two posi- 
tions symmetrically placed on either side of the central 
star. In this case the slit was 2 arcsec wide by 10 arcsec 
long and was centered 3 arcsec to the E and then 3 arcsec 
to the W of the central star with a slit position angle of 
90°, as shown in Figure [6] As with Tc 1 these slit posi- 
tions are well inside the outer regions of the nebula. The 
combined exposure times at each slit position for each 
grating setting were of order 60 min. We flux calibrated 
the NGC 6543 spectra and then averaged together the 
two slit positions to get a final spectrum interpolated for 
the line of sight to the central star in the same way as 
was done with the LCO spectra. We determine the ab- 
solute spectrophotometry of our calibration to have an 
accuracy of better than 7 percent. In Figures [7][9] we 
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TABLE 4 

Optical Forbidden Emission Line Fluxes 



Species 


A(A) 


U P 9 1 QQ 

rrez- loo 




NGC 6543 


Tc 1 




F c 


F 


F c 


F 


F c 


HI 


4861 


3.49(-12) 


1.26(-11) 


1.51(-11) 


2.10(-11) 


1.58(-12) 


3.37(-12) 


C I 


4622 


2.76(-16): 


1.07(-15): 










C I 


8727 


1.63(-15) 


2.86(-15) 






<1.97(-16) 


<2.76(-16) 


C I 


9824 


1.16(-15): 


1.84(-15): 










C I 


9850 










<9.74(-16) 


<1.28(-15) 


N I 


5198 


3.08(-14) 


9.97(-14) 


4.37(-15): 


5.92(-15): 


3.79(-16): 


7.61(-16): 


N I 


5200 


1.88(-14) 


6.09(-14) 


3.92(-15): 


5.30(-15): 


4.04(-16): 


8.10(-16): 


N II 


5755 


3.60(-14) 


1.01(-13) 


4.70(-14) 


6.13(-14) 


2.00(-14) 


3.68(-14) 


N II 


6548 


3.08(-12) 


7.37(-12) 


7.66(-13) 


9.59(-13) 


6.17(-13) 


1.04(-12) 


Nil 


6583 


9.38(-12) 


2.23(-ll) 


2.55(-12) 


3.18(-12) 


1.92(-12) 


3.21(-12) 


OI 


5577 


6.33(-16) 


1.84(-15) 






6.61(-17): 


1.25(-16): 


I 


6300 


7.36(-14) 


1.85(-13) 


2.12(-14) 


2.69(-14) 


2.45(-15) 


4.24(-15) 


O I 


6364 


2.66(-14) 


6.60(-14) 


7.03(-15) 


8.82(-15) 


9.68(-16) 


1.66(-15) 


O II 


3726 


4.51(-13) 


2.25(-12) 


1.40(-12) 


2.12(-12) 


1.68(-12) 


4.36(-12) 


O II 


3729 


2.00(-13) 


9.97(-13) 


6.41(-13) 


9.70(-13) 


1.11(-12) 


2.90(-12) 


O II 


7320 


5.02(-12) 


1.06(-13) 


1.74(-13) 


2.10(-13) 


1.21(-13) 


1.89(-13) 


O II 


7330 


4.74(-14) 


9.97(-14) 


1.79(-13) 


2.17(-13) 


1.02(-13) 


1.59(-13) 


OIII 


4363 


<6.29(-16) 


<2.64(-15) 


3.73(-13) 


5.40(-13) 


7.97(-15) 


1.87(-14) 


III 


4959 


5.72(-16) 


1.99(-15) 


4.15(-11) 


5.73(-ll) 


6.72(-13) 


1.41(-12) 


O III 


5007 


2.72(-15) 


9.34(-15) 


1.24(-10) 


1.71(-10) 


2.00(-12) 


4.17(-12) 


p II 


4669 


<6.29(-16) 


<2.40(-15) 






<2.07(-16) 


<4.58(-16) 


p II 


7876 


1.20(-16) 


2.34(-16) 






<5.98(-17) 


<8.88(-17) 


S I 


4589 


<8.88(-16) 


<3.48(-15) 






<2.07(-16) 


<4.66(-16) 


S I 


7725 










<7.99(-17) 


<1.20(-16) 


S II 


4069 


1.23(-13) 


5.64(-13) 


6.84(-14) 


1.01(-13) 


8.36(-15) 


2.06(-14) 


S II 


4076 


4.47(_14) 


2.04(-13) 


1.24(-14) 


1.83(-14) 


2.40(-15) 


5.91(-15) 


S II 


6716 


4.69(-13) 


1.09(-12) 


1.38(-13) 


1.71(-13) 


4.50(-14) 


7.43(-14) 


SII 


6731 


9.77(-13) 


2.26(-12) 


2.54(-13) 


3.15(-13) 


7.15(-14) 


1.18(-13) 


S III 


6312 


2.24(-15) 


5.60(-15) 


1.71(-13) 


2.17(-13) 


9.22(-15) 


1.59(-14) 


S III 


9069 


3.44(-13) 


5.85(-13) 


4.91(-12) 


5.63(-12) 


3.19(-13) 


4.38(-13) 


CI III 


5518 


4.57(-16) 


1.35(-15) 


6.54(-14) 


8.65(-14) 


5.06(-15) 


9.62(-15) 


CI III 


5538 


7.90(-16) 


2.32(-15) 


9.27(-14) 


1.22(-13) 


5.39(-15) 


1.02(-14) 


Ar III 


5192 






1.05 (-14) 


1.42(-14) 


5.29(-16) 


1.06(-15) 


Ar III 


7136 


2.78(-15) 


6.02(-15) 


3.56(-12) 


4.35(-12) 


1.53(-13) 


2.42(-13) 


Ar III 


7751 


3.75(-16) 


7.34(-16) 


8.30(-13) 


9.89(-13) 


3.85(-14) 


5.78(-14) 


Ar IV 


4711 






1.77(-13) 


2.49(-13) 






Ar IV 


4740 






2.08(-13) 


2.92(-13) 






Fell 


4244 


<1.44(-15) 


<6.28(-15) 


1.10(-15): 


1.61(-15): 


<3.12(-16) 


< 7.47(-16) 


Fell 


4359 


3.56(-15) 


1.50(-14): 


1.39(-15): 


2.01(-15): 


<2.08(-16) 


<4.87(-16) 


Fell 


5159 


1.58(-15) 


5.17(-15): 










Fell 


8617 


2.75(-15) 


4.91(-15): 






<1.98(-16) 


<2.78(-16) 


Ni II 


6667 


1.49(-15): 


3.48(-15): 


4.69(-16): 


5.83(-16): 


2.29(-16): 


3.80(-16): 


Ni II 


7376 


1.71(-15): 


3.58(-15): 






<1.00(-16) 


<1.55(-16) 



a In units of erg s 1 cm 2 , numbers in parentheses are exponents, colons after values indicates uncertain detections. 

b Observed flux. 

c Extinction-corrected flux. 



show portions of the nebular spectra of the three PNe 
that were used in the analysis of the emission lines and 
which show both the strong diagnostic forbidden lines 
and the weaker recombination lines of O II. 



4.2. Nebular Emission Line Intensities 

The emission line fluxes have been measured from the 
final co-added and averaged spectra using the IRAF splot 
routine. The measurements were straightforward be- 
cause few of the lines of interest showed evidence of sig- 
nificant blending. The resultant observed intensities for 
He2-138, NGC 6543, and Tc 1 are given in Table H for 
lines that can be used to obtain T e , n e , and extinction 
along the lines of sight. The observed intensities have 
been corrected for extinction by taking the flux ratios of 
multiple unblended Balmer and Paschcn line pairs from 
the same upper levels and determining the logarithmic 



extinction at H/3, Cnp, from the expression 

c H p = [X H p/{Xi - X 2 )] x log 10 {A 1 F 2 X 2 /A 2 F 1 X 1 ) , 

(3) 

where A\ >2l A1.2, and F\ t2 are the spontaneous emission 
coefficients, wavelengths, and observed fluxes for a spe- 
cific Balmer and Paschen line pair, and X±_2.Hd a re the 
galactic extinction law values fitted bv iHowarthl (|1983l ) 
at the wavelengths of the lines and H/3 respectively (as- 
suming R — 3.2). Individual emission line fluxes were 
then corrected using the relation 

F C (A) = 10 C ^ X < A )/ A '^F(A) , (4) 

where F c is the corrected flux. Taking the average of val- 
ues obtained from multiple line pairs for our lines of sight 
we derive values of Ch/3 = 0.56, 0.14, and 0.33 for He2- 
138, NGC 6543, and Tc 1 . These values are in good 
agreem ent with those of ICahn. Kaler. fc St anghelliui 
(|1992f) . who obtained global values of 0.40, 0.12, and 0.28 
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Fig. 10. — Diagnostic diagram for He2-138. Dashed lines indicate 
n e diagnostic curves and solid lines T e diagnostics curves derived 
from emission line intensities listed in Table [4] input into the ratios 
listed in Table E] 
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Fig. 11. — Diagnostic diagram for Tc 1. Dashed lines indicate 
n e diagnostic curves and solid lines T e diagnostics curves derived 
from corresponding emission line intensities listed in Table [4] input 
into the ratios listed in Table [6] 

for the three PNe. The extinction corrected fluxes, F c , 
are listed in column 4 of Table [H including the upper lim- 
its to fluxes of undetected lines, which have been taken 
to be the 3cr rms flux of the noise of the neighboring 
continuum. 

4.3. Plasma Diagnostics and Emission Measures 

Emission measures and relative abundances of ions are 
normally determined from their forbidden line intensi- 
ties, which have a sensitive dependence upon the ki- 
netic temperature and density of the gas. Electron tem- 
peratures T e have been determined for regions of dif- 
ferent ionization primarily from the ratio of auroral to 
nebular line intensities of [O I], [S II], [N II], [O II], 
[O III], [S III], and [Ar III] using complete radiative 
and collisional multi-level cal culations similar to th ose 
in the IRAF nebular packag e dghaw fc Dufounll995h . as 
described by ISharpee et al.1 (|2007l ) in their study of s- 
process elements in PNe. Similarly, electron densities n e 
are sensitive to certain line ratios such as [O II] A3726/29, 
[S II] A6716/31, [CI III] A5518/38, and [Ar IV] A4711/40. 
We have used all of these forbidden line ratios, when they 
were observed, and corresponding atomic data listed in 
Table [5] to calculate appropriate values of T e and n e in 
Figures [T0HT21 The resulting values of T e and n e are 
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Fig. 12. — Diagnostic diagram for NGC 6543. Dashed lines in- 
dicate n e diagnostic curves and solid lines T e diagnostics curves 
derived from corresponding emission line intensities listed in Ta- 
ble [4] input into the ratios listed in Table [6] 

listed in Table [5] with their formal uncertainties. The 
densities and temperatures derived from the different 
lines are generally consistent with each other with the 
exception of the [N I] for NGC 6543 and Tel, and the 
[S II] density for He2-138, as is evident from the plots 
in Figures [T0rfT2l The disparate densities deduced for 
He2-138 from the different lines could be real; the result 
of inhomogeneities. The [N I] lines, on the other hand, 
are very weak and thus the densities from that doublet 
are very uncertain. 

In order to properly account for all relevant physi- 
cal processes when converting the observed emission line 
flux F c into the emission measure we consider here the 
full definition of the emission measure. The extinction- 
corrected flux of an optically thin emission line along a 
line of sight is 

F c = e 2 [hv /(Air)] J n u A ul dl , (5) 

where 9 2 is the angular area of the gas being observed, 
and A u i and n u are the line transition probability and 
number density of the upper level. The stronger forbid- 
den transitions normally have direct collisional excitation 
from the ground state as the predominant mechanism 
exciting the line, therefore it is convenient to write the 
equation of statistical equilibrium governing the popula- 
tion of the upper level in terms of the ion and electron 
densities as 

uk — ^e^lQlu {T e )[l + Z(n e ,T B ,J v )\, (6) 

where qi u (T e ) is the collision coefficient between the 
ground state and upper level, and the term £(n e , T e ,J v ) 
represents all other processes contributing to the popu- 
lation of the upper level, e.g., radiative cascading and 
collisional population from upper levels, collisional de- 
excitation to lower levels, continuum and resonance fluo- 
rescence, and recombination. We formulate the equation 
this way in order to isolate the term n e ni, whose integral 
along the line of sight is the emission measure of the ion. 
For the processes listed above, the expression for £ is 

I k>u k>l.k^u 
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TABLE 5 

Atomic Data References 



Species 



Transition Probabilities 



Collision Strengths 



C 1 

N I 
N II 

O I 

O II 
O III 

p II 

S II 
S III 



CI III 
Ar III 

Ar IV 



Fe II 

(159 levels) 



Ni II 

(76 levels) 



Nussbaumer & Rusca (1979) 
Froese Fischer fc Saha (19851 

Zeippen (19821 

Froese Fischer fc Tachiev ( 20041 
^^^^^^^^^^^^^^a'^^TO^ - 

Wiese. Fuhr. fc Deters 119961 



Deters 



Baluia & Zeippen (1988) 
Mendoza (19831 

Zeippen (19821 
Wiese, Fuhr, fc Deters f 19961 
Nussbaumer & Storey (1981) 
Wicse^F^uhr^ " 



Kaufman fc Sugar (19861 
Mendoza fc ZeipperT7 l982b1 
Mendoza fc Zeippen (1982a) 
j<raman_et_al ii (T993) 
\/ejmer i ^ejmer J _fc Fcrland (1996) 
Mendoza fc Zeippen (1982bT ~ 
Fiejse J _Sjnitli J _&_Cala^^^S95 l ) 
LaJohn k Luke (1993^ 



Mendoza & Zeippen (1982a) 
Mendoza fc Zeippen (19831 

Mendoza fc Zeippen (1982a) 
Kaufman fc Sugar ( 19861 
^Jus^bairmCT_&Stoj^^"Q988) 
Garstang (19621 
Nahar (19951 

^clmabeL^^Schultz-Johanning, 
Bautista fc KallmanT 20011 
Nussbaumer fc Storey (19821 
Kurucz (1992) 



Kock (2004 



Pcauignot & Aldrovandi (1976) 
Thomas fc Nesbet (19751 
i Jol i mson Ji Fjurke J _& Kingston (1987) 
Berrington fc Burke (19811 
Dopita. Mason, fc Robb (19761 
Stafford et al. (19941 
Saraph, Scaton, & Shcmming (1969) 
Lennon fc Burke (19941 
Berrington & Burke (1981) 
Berrington (1988) 
LeDourneuf & Nesbet (19761 
Pradhan (19761 
McLaughlin fc Bell (19931 
Aggarwal (19831 

Aggarwal, Baluia. fc Tullv (19821 
Balui a. Burke, fc Kingston (1980) 
Baluia, Burke, & Kingston (1981) 
Lennon fc Burke (19941 
Taval (20041 
Kruger fc Czvzak (1970 1 
Kccnan ct al. (1996) 
Mendoza (T9821 



Butler fc Zeippen (19891 
Johnson & Kingston (1990) 
Kruger fc Czvzak (197pT ~ 
^ejrJP^n^IjeFtorlotj^ Butler 




Bautista (2004 



TABLE 6 

Electron Temperatures and Densities 



Diagnostic 



Hc2-138 



NGC 6543 



Tc 1 



[N I] A5198/A5200 
[S II] A6716/A6731 
[O II] A3726/A3729 
[CI III] A5518/A5538 
[Ar IV] A4711/A4740 

[O I] (A6300+A6364)/A5577 

[S II] (A6716+A6731)/(A4069+A4076) 

[O II] (A3726+A3729)/(A7320+A7330) 

[N II] (A6548+A6583)/A5755 

[S III] A9069/A6312 

[Ar III] (A7136+A7751)/A5192 

[OIII] (A4959+A5007)/A4363 



Density (cm 3 ) 
7000t^ 000 
15000 

+9000 



7500 
7500 H 



3000 



Temperature (K) 
hl200 



9000Z 



7000 
6500 
6000 



6000 

+3000 



J i 

+ 000 
700 
+500 
300 



goQ+900 
yuu -500 
5000+ 17000 

JUUU -3000 

eooot^oo 

5000S 
4500™ 



9000 
12000 
10300 
8500 
8000 
8200 



7000 
5000 
+ 3000 
2000 
+ 800 
700 
+ 500 
400 
+400 
300 
+ 200 
200 



400600 

+ 3800 
-1200 
+ 800 
-000 
1MIO 



3000 
2000 
3000 H 



140001?™ 

-,+5000 
-3000 
-,+2700 
-1700 
t-600 
-500 

95001™ 

h800 
-500 
-,+500 
-400 



90001 
10500j 
8500^ 
9500^ 
9000;! 
9000j 
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^ t Quk + y t n kJkuBku 
-n u JukB u k + n e n l+1 a u > . (7) 

k^u J 

Here, Jku and B^u are the mean intensity and Einstein 
B coefficient for radiative (de-)excitation from level k to 
level u, and a u is the effective recombination coefficient 
into level u. 

For most strong forbidden lines direct excitation from 
the ground states predominates, and £ < C 1. How- 
ever, for certain ions, e.g., Fe II a nd Ni II (iLucvl 119951; 
Bautista. Peng &; Pradhanl Il996t iBautista fe Pradhanl 
19961: i Bautistal l2004f ). and certain transitions of CNO 
(jGrandil Il975[ ) other processes such as radiative exci- 
tation and strong coupling between other excited levels 
contribute to some of the stronger forbidden transitions, 
such that £ > 1 for these lines under certain conditions. 
The intensities of these lines do not retain a simple linear 
dependence on n e n\ , and it is important to treat their ex- 
citation via detailed multi-level calculations that involve 
the radiation field. 

We have used detailed calculations of the relevant level 
populations for all of the ions using the values of T e , n e , 
and J v appropriate for the level of ionization and aug- 
mented by incorporation of additional levels and pro- 
cesses for the ions, to determine the emission coefficients 
Qiu(Te) and £(rt e , T e , J„) for all of the lines listed in Ta- 
ble HI The radiation fields J v for our slit positions have 
been taken from observations of the central stars by IUE 
and FUSE for frequencies below the Lyman limit. For 
frequencies above the Lyman limit we have extrapolated 
the observed stellar continua by assuming a black body 
flux at the appropriate temperature for the central stars. 
The dilution factors at our slit positions were rather small 
so that radiative excitation by the stellar continua was 
not competitive in the population of any level that we 
considered. For Fe II and Ni II we have used the ex- 
plicit multi-level population processes and cross sections 
of Bautista, Pradhan, and collaborators to calculate line 
strengths for these ions. All known processes that might 
make significant contributions to the line intensities have 
been included in the above calculations, and correspond- 
ing values of T e and n e have been used that are appro- 
priate for the ionization state for each line. Reddening 
corrected fluxes have then been used to compute ground 
state emission measures EMi for those ions for which UV 
absorption was also observed. 

Combining equations [50 produces the following ex- 
pression for the emission measure of an ion i, 

EMi = J n e nidt 

= AttF c ^ A uk/ {hv A u i8 2 q lu (T e ) 

k<u 

x [l+£(n e ,T e , J„)]} . (8) 

The extinction-corrected fluxes from Tabled! the atomic 
data and cross sections from the references in Table [5l 
and the temperatures and densities listed in Table [6] ap- 
propriate to the different ions depending upon their level 
of ionization have all been used to determine the emission 
measures of ions from equation [8] using the intensities of 



TABLE 7 
Forbidden Line Emission Measures 



Species 


A(A) 


EM, 






(cm -6 pc) 




Hc2-138 


CI 


8727 


0.404 


PII 


7875 


<0.025 


Fell 


5159 


0.61 


Ni II 


7378 


0.0124 


OI 


6300 


12.36 




6364 


13.79 


SII 


4069 


13.99 




4076 


14.92 




6716 


11.67 




6731 


12.84 


sin 


6312 


3.93 




9069 


3.89 


NGO 6543 


NI 


5198 


0.029 




5200 


0.031 


I 


6300 


0.227 




6364 


0.234 


SII 


6716 


0.202 




6731 


0.201 


S III 


6312 


5.493 




9069 


5.598 


Tc 1 


OI 


6300 


0.060 




6364 


0.073 


SII 


4069 


0.099 




4076 


0.084 




6716 


0.136 




6731 


0.139 


S III 


6312 


0.736 




9069 


0.990 



the various lines for our sample of PNe. The resulting 
emission measures are presented in Table [7] 

5. COMPARATIVE COLUMN DENSITIES FROM 
EMISSION AND ABSORPTION 

5.1. Correction for Different Lines of Sight 

In order to compare directly the results of the absorp- 
tion and emission abundance analyses the emission mea- 
sures have to be converted to effective column densities, 
or vice versa, by dividing the emission measures by the 
electron density appropriate for each ion. If we designate 
<n e >i as the mean value of n e in the emitting region of 




Fig. 13. — A schematic indicating the relevant angles discussed 
in the text that apply the derivations of the volume and shell cor- 
rection factors £ vo i and Cshell, assuming that the nebula is a perfect 
sphere. 
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the ion i, the equivalent emission column density for the 
ion can be written as 



N em = EMi/(( <n e >i) 



(9) 



since the integrals that define the emission measure and 
the column density of an ion differ only by the factor 
of the electron density in the emission measure. The 
constant C is a normalization factor that corrects for the 
different path lengths along the two lines of sight, viz., 
our emission line of sight passes through the entire nebula 
whereas the absorption spectrum line of sight penetrates 
only to the central star. A derivation of C is given below 
that allows for the fact that our observations record the 
flux within a rectangle subtended by the entrance slit of 
the spectrograph, and this slit is offset from the center 
of the nebula to avoid contamination of the spectrum by 
the central star. 

A generic representation of our emission-line measure- 
ment geometry is shown in Figure I13[ with the ideal- 
ization that the appearance of the nebula is perfectly 
round in the sky. We consider two fundamentally differ- 
ent simplifications for the distribution of any particular 
ion within the nebula. The first such representation is 
a uniformly filled sphere. If the height of the slit y is 
longer than the chord through the nebula, the volume 
element V of the sphere that is interior to a projection 
of the slit with a width equal to X = 2 — 0\ is given by 



(10) 



If we consider that the integration range in equation [5] 
is over a distance equivalent to the angle subtended by 
r , we can define the correction factor Cvol,o to be sim- 
ply V derived above normalized to the volume within a 
rectangular box of dimensions x y r ; hence 

TT [Of - 0\ + 30 2 r (0 2 - 00] 



c 



V 



vol.O 



30 x y 



(11) 



If y does not subtend the entire chord, we must reduce 
Cvoi.o by a factor 



Fvol = - 
7T 



_» 4 «. 



2 sin 



28; 



where 



(12) 



(13) 



is the radius of the circle that represents the intersection 
of the surface of the sphere with a plane that is aligned 
with the center line of the slit. The factor F vo i is based on 
the approximation that 8 X <C Z , since it represents the 
area subtended by lines bounded by y inside the circle 
relative to the total area of the circle, but only for a circle 
coincident with the slit centerlinc. The final value for Cvol 
which applies to a fully filled sphere is given by 

Cvol = Fvol Cvol.O (14) 

Our second representation differs from the first in that 
the material is assumed to be distributed in a thin shell, 
rather than throughout the entire volume of the nebula. 
In a development similar to the one that we performed 



TABLE 8 
Angles and Correction Factors 

Nebular Identification 
Quantity a He2-138 NGC 6543 Tel 





3.5 


9.8 


4.8 




2.0 


2.0 


1.9 


e-> 


3.5 b 


1.0 


3.9 


By 


4.0 


10. 


4.0 


Cvol 


0.985 


1.80 


1.48 


Cshell 


4.12 


2.22 


2.64 



Angles given in arc seconds. 
One side of the slit extended beyond the edge 
of the nebula (#2 — 4.0), hence 62 is set to 8 r . 

for the volume-filled nebula, we compute the area A of a 
projection of the slit on the surface of the sphere, under 
the condition that it is long enough to cover the entire 
chord, 

A = 2tk0 t (O2 — Ox) ■ (15) 

In this case, the normalization box has the dimensions 
x y multiplied by the thickness of the shell. Since the 
appropriate path length for equation [5] in this case is the 
shell thickness, this thickness cancels out in the equation 
for shell correction factor Csheii,o 5 leaving us with the ex- 
pression 

A 2-K0 r {0 2 -0 1 ) 



c 



shell, — 



(16) 



The reduction factor -F s heii for y < 20 z is given by 
2 . 



K 



hell 



20 2 



(17) 



and this factor reverts to unity for y > 20 z . As before, 

Cshcll = -Fshcll Cshcll.O (18) 

Values for the angles and correction factors Cvol and 
Cshcll are given in Table [5J and are based upon nebular 
diameters taken from the literature, as discussed below. 
Both values of C for NGC 6543 and Tc 1 are not very far 
from 2 because the slit heights were only about half the 
nebular diameters and the slits were positioned rather 
close to the central stars. The relative change in going 
from Cvol to Csheii is large for He2-138 because the slit 
size was comparable to the nebular diameter, and it was 
positioned near the nebula's edge. 

The images of He2-138, NGC 6543, and Tc 1 in Fig- 
ure[5]show that all three PNe possess a degree of spherical 
symmetry for the overall structure, but also have embed- 
ded asymmetrical, inhomogeneous features, e.g., clumps 
and possible bipolar structure. Differences between the 
emission and absorption lines of sight therefore depend 
not just on the footprint of the spectrograph slit on the 
nebulae and whether the nebulae can be represented as 
filled shells or thin rings, but also upon small-scale inho- 
mogeneities that lie along one line of sight but not the 
other. Small-scale structure could well be the dominant 
cause of differences between the separate lines of sight, 
and such structure commonly depends on the level of 
ionization. 

Our emission spectra provide information on the ex- 
tent of both large-scale geometrical and small-scale in- 
homogeneity effects from the intensity variations of the 



Independent Emission and Absorption Abundances 



15 



TABLE 9 

Comparitive Column Densities from Emission and Absorption Lines 



Species 


log N abs 




< n e >i 


log N em a 


log N abs - log N em a 




(cm -2 ) 


(cm -6 pc) 


(cm -3 ) 


(cm -2 ) 










He2-138 (C 


ml — 0.985 Cshcll 


= 4.12) 




Cl 


13.93±0.16 


0.404 


7000 


(14.26, 13.64) ±0.31 


(-0.33, 0.29) ±0.35 


PII 


13.58±0.21 


<0.025 


7000 


<(13.05, 12.43) ±0.42 


>(0.53, 1.15) ±0.47 


Fell 


14.37±0.14 


0.61 


7000 


(14.44, 13.82) ±0.64 


(-0.07, 0.55) ±0.66 


Ni II 


12.77±0.17 


0.0124 


7000 


(12.74, 12.12) ±0.21 


(0.03, 0.65) ±0.27 


OI 


>15.48±0.37 


12.8 


7000 


(15.76, 15.14) ±0.34 


>(-0.28, 0.34) ±0.50 


SII 


>15.49±0.09 


12.3 


10000 


(15.58, 14.96) ±0.26 


>(-0.09, 0.53) ±0.28 


S III 


15.11±0.13 


3.89 


7500 


(15.21, 14.59) ±0.11 


(-0.10, 0.52) ±0.17 






NGC 6543 (Cvoi = 1-80 Q hcl 


1 = 2.22) 




Ol 


13.48±0.51 


0.23 


5000 


(13.89, 13.80) ±0.25 


(-0.41, -0.32) ±0.57 


NI 


12.29±0.25 


0.030 


4000 


(13.10, 13.01) ±0.21 


(-0.81, -0.72) ±0.33 


SII 


13.64±0.07 


0.20 


5000 


(13.83, 13.74) ±0.14 


(-0.19, -0.10) ±0.16 


S III 


15.22±0.14 


5.56 


5000 


(15.27, 15.18) ±0.04 


(-0.05, 0.04) ±0.15 






Tc 1 (Cvoi = 1-48 (shell = 


2.64) 




Ol 


13.67±0.06 


0.064 


2000 


(13.82, 13.57)±0.26 


(-0.15, 0.10) ±0.27 


SII 


13.67±0.09 


0.137 


3000 


(13.98, 13.73)±0.14 


(-0.31, -0.06) ±0.17 


sin 


14.93±0.07 


0.92 


3000 


(14.80, 14.55) ±0.04 


(0.13, 0.38) ±0.08 



a Numbers in parentheses show the outcomes for £ vo i and Caheii, respectively. These are followed by the error limits 
arising from measurement uncertainties. 



emission lines along the slit length. The slit lengths we 
used were 4 and 10 arcsec with a spatial resolution of 
1 arcsec along the slit, so we have relatively few inde- 
pendent spatial resolution elements. Nevertheless, the 
range in distance of the slit from the central star over its 
length is comparable to the offset of the slit center from 
the central star. Thus, variations of intensity along the 
slit due to the overall geometry of the nebulae should be 
comparable to the intensity differences due to the differ- 
ent path lengths of the absorption and emission sight- 
lines. We have measured the variations in intensity of 
the [S III] and [O I] lines, representing our highest and 
lowest ionization species, along the slit in the three PNe. 
We find for Tc 1 that both the [S III] and [O I] lines 
have a very uniform distribution of intensity throughout 
the full slit, and with no measurable differences between 
the two slit positions. Thus, for Tc 1 the measurements 
indicate that the emission and absorption lines of sight 
are likely to be very similar. 

For He2-138 the [S III] and [O I] lines have virtually 
identical, smooth intensity distributions where the inten- 
sity peaks at the center and decreases outward toward the 
ends of the slit where it falls off rapidly near the ends. 
The intensity profile is more characteristic of a filled vol- 
ume than a thin shell distribution of gas, but the smaller 
size of this nebula causes it to fill only 3 arcsec of the 
4 arcsec slit length so the geometrical factor zeta repre- 
sents an important correction for this object. Since there 
is no indication of differences in the spatial distribution 
of the ions based on their ionization level the geometrical 
normalization factor zeta is the same for all of the lines. 
The fact that the emission spectra sampled the outer 
edge of the nebula makes the corresponding geometrical 
correction for this PN rather uncertain, as was explained 
in §4.1.11 Thus, we consider the results for He2-138 to 
be less reliable than those of Tc 1 and NGC 6543. 

NGC 6543 presents a more complicated picture in 
terms of the differences between the two lines of sight. 
The [S III] completely fills the slits with a uniform inten- 
sity for one of the slit positions, but shows variations of 



^25 percent in the other position. The [O I] completely 
fills the slits also, but shows large variations along the slit 
near the center. Thus, the lowest ionization species in 
this PN display a pronounced small-scale structure that 
may cause the emission and absorption lines of sight to 
be quite different for the lowest ionization lines. Based 
on the intensity variations one must admit the possibility 
of differences in the column densities along different lines 
of sight for the neutral species to be as large as factors 
of 3 for this object - an uncertainty that compromises its 
usefulness for the neutral species [O I] and [N I] . 

If the stellar absorption and nebular emission spec- 
troscopy are obtained at sufficiently high resolution one 
can use radial velocity information from resolved line pro- 
files to match velocity components of optical emission 
with the corresponding UV absorption produced in the 
same velocity intervals. A comparison of these quanti- 
ties within the same velocity interval of the gas provides a 
more accurate assessment of the comparative abundances 
than comparing the total emission measure and column 
density integrated over the full profiles. Local values of 
the emission column density can be determined for spe- 
cific kinematic regions within the nebulae. Averaging 
these values over the full velocity range for the nebu- 
lar shell will produce the global emission column density 
for the ion, and also information on its fluctuations as 
a function of velocity. Since thermal and expansion ve- 
locities of nebulae are of order 10-20 km s _1 a spectral 
resolution less than 10 km s _1 is optimal to perform the 
analysis this way. Our emission spectra lack the neces- 
sary spectral resolution to perform such an analysis, and 
therefore we work with the integrated (over wavelength) 
emission measures and column densities. 

5.2. Comparison of Absorption & Emission Column 
Densities 

The absorption column densities obtained from UV 
resonance lines refer specifically to those ions occupying 
the lower level of the transition. In order to obtain the 
total column density of the ion the column densities for 
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all the individual fine-structure levels of the ground state 
must be summed together. In our STIS spectra some 
ions had blended absorption profiles for one or more of 
the transitions from the ground state fine-structure levels 
that prevented us from deriving the column densities for 
those levels. We have determined the column densities of 
ions in those levels by taking values of T e and n e obtained 
from the emission lines for that ion to solve for the level 
populations relative to the levels for which column den- 
sities were determined. Additionally, when more than 
one resonance multiplet of an ion has yielded a column 
density we have computed the mean value for the ion 
by weighting individual values according to the inverse 
square of their uncertainties. These calculations, which 
have been applied to the absorption column densities in 
Table [H yield the total ion column densities, N a b s , to 
the central star, and these are listed in Table ED together 
with the formal errors that result from the quantifiable 
uncertainties that are discussed below. 

Emission measures for the same ions that have been ob- 
served in absorption, and which appear in Table[7]for our 
sample of PNe, are also presented in Table E3 For ions 
where more than one forbidden line yields an emission 
measure we have determined the average of the values, 
with stronger weight being given to the lines of higher 
intensity and lower Boltzmann factor. Using the corre- 
sponding values of the electron density for each of the 
ions the resulting emission column densities, N em , have 
been determined from the emission measures from equa- 
tion [9] These are given in the penultimate column of 
Table E3 together with the combined errors, having been 
normalized to the absorption path length by dividing by 
the geometrical factor £. 

A comparison of the values of N a t, s and N em for the 
different ions from the two completely independent abun- 
dance analyses shows moderately good agreement, with 
the exceptions of P II in He2-138 and N I in NGC 6543. 
Absolute abundances determined from the forbidden 
emission lines and UV absorption lines give the same 
results within ±0.3 dex for adjacent lines of sight, which 
is comparable to the combined formal errors of the anal- 
yses. The la errors in the column densities derived from 
the analyses represent the uncertainties that are quan- 
tifiable. There are also systematic errors that arise from 
assumptions rather than measurement uncertainties, and 
both sources of error affect the accuracy of our compari- 
son of forbidden and recombination line column densities. 

The primary sources of error for the absorption col- 
umn densities are (a) the determination of the proper 
continuum level, (b) the low S/N of the intensities of 
weak absorption lines, (c) the insensitivity of intensity 
to column density for saturated lines, and (d) the deter- 
mination of total ground state column density for states 
with fine-structure levels when absorption from one or 
more of the levels is either not observed or saturated. 
The main sources of error for the emission column den- 
sities are uncertainties in (a) the flux calibration of the 
echelle spectra, which are at the 5-10 percent level, (b) 
collision strengths for some of the forbidden lines, and (c) 
the correct values of T e and n e that correspond to each of 
the transitions, as assigned to the various ions from the 
diagnostics shown in Figures [IT)lTl"2l The atomic data 
for most of the forbidden lines that we have used for di- 
agnostics and the determination of column densities are 



believed to be known to better than 30 percent accu- 
racy. With the exception of the [P II] line the current 
values for most of the forbidden line collision strengths 
and transition probabilities are the result of calculations 
by independent methods over the past three decades that 
have converged on values that are in good agreement with 
each other and that have changed little over the past five 
years. Thus, the atomic data are not likely to be major 
sources of error. Rather, the largest sources of formal 
errors in the emission column densities are uncertainties 
in the values of temperature and density. Because line 
intensities depend upon these two parameters, errors in 
T e and n e translate to errors in the column density. Most 
of the lines are in the low density limit and therefore the 
emission measures are rather insensitive to density. How- 
ever, because of the Boltzmann factor the line intensities 
are sensitive to T e . The errors caused by uncertainties 
in the temperature, together with uncertainties in inten- 
sity measurement and flux calibration, form the basis for 
the combined error that is presented in Table for each 
emission column density. To these uncertainties must 
be added the unknown errors in collision strengths and 
those differences that small-scale inhomogeneities may 
cause between the lines of sight. 

Several features of Table [9] merit comment. First, due 
to a combination of weak, saturated, or strongly blended 
UV absorption lines coupled with the failure of our nebu- 
lar spectra to detect forbidden lines from some ions, there 
are relatively few ions for which we were able to derive 
independent abundances from both UV absorption and 
forbidden emission lines. Even with the relatively long 
slit used to sample substantial portions of the PNe shells 
we were not successful in detecting weak emission lines 
from a number of the ions for which column densities had 
been measured from the STIS spectra. 

Second, with the exception of S +2 all of the ion species 
listed in Tables [7] and [9] are the lowest ionization stages 
that have ionization potentials greater than 13.6 eV. This 
means that some fraction of most of the ions we have 
measured could exist in cold, neutral gas residing either 
within dense clumps embedded inside the nebula or in 
a foreground shell of material around the nebula. Such 
gas would complicate our analysis by increasing the ab- 
sorption column density without having any effect on 
the emission lines, leading to legitimate differences be- 
tween the emission and absorption column densities. For- 
tunately, we can test for this possibility by comparing 
N a bs(0 I) with jV a ;, s (S II). The ionization fraction of O 
is closely coupled to that of H through a charge exchange 
react i on that has a large rate constant (|Field fc Steigmanl 
ll971tlChambaud et alJ ll980). which guarantees that the 
amount of O I in the ionized gas is quite low and that 
practically all of the O is neutral in H I gas. By contrast, 
in an H II region a reasonable fraction of the S will be 
in the form of S + since its ionization potential is high 
(23.3 eV). Furthermore, in an H I region nearly all of the 
S should also be singly ionized. Therefore, from H I gas 
we expect to find the ratio N a b s (0 l)/iV a & s (S II) to be 
approximately equal to the solar value of [O/S] = 1.46, 
assuming that neither of the elements are significantly 
condensed onto dust grains nor are they enriched or de- 
pleted by nuclear processes within the AGB progenitor of 
the central star. A ratio smaller than this value signifies 
progressively less contribution to the column densities 
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from neutral gas. 

For Tc 1 we have found that N abs (0 I)«A a6s (S II), 
which indicates that any contribution from neutral ma- 
terial must be so small that it can be neglected for 
our study. The situation for NGC 6543 is not quite so 
straightforward because our inferred value of N abs (0 I) 
for the nebula is based on the marginal detection of O I*. 
Our ability to directly measure O I in the ground fine- 
structure level is compromised by possible P II A1301.87 
absorption from foreground ISM gas at a velocity v = 
-14 km s _1 , which appears at the same wavelength as 
the velocity-shifted O I A1302.17 line from the nebu- 
lar shell. Nevertheless, we can derive an upper limit to 
the O I column density for the nebula from this fea- 
ture, which has equivalent width EW — 10 mA, by as- 
suming the foreground P II absorption to be negligible. 
When we do this, we derive a value \ogN(0 I) = 13.2, 
an amount that is above the lower bound for our calcu- 
lated log N abs (0 I) that is listed in Table [9] However, 
this value is still substantially lower than our measure- 
ment of N abs (S II), so once again we are assured that 
contamination of the column densities from neutral gas 
is negligible for NGC 6543. 

We are unable to make any assertion about 
N a bs(0 I)/N abs (S II) toward He2-138 because both col- 
umn densities were recorded as lower limits (the lines are 
strongly saturated; see footnote c in Table [3J). However, 
in the spectrum of the central star for this nebula we see 
absorption features from excited H2 at v = -62 km s" 1 , 
which is a strong indication that we are viewing a pho- 
todissociation region at the inner edge of a neutral shell 
surrounding the nebula. Thus, it is possible for this one 
object that N abs for ions that can exist within H I regions 
could add to the contributions from the ionized nebula. 
This may explain why N abs (P II)»7V em (P II) for this 
PN, although it is then puzzling why the discrepancies 
for Fe II and Ni II are not nearly as large unless both of 
them are condensed onto grains. 

Finally, we point out that heavy element recombina- 
tion lines for the ion species that we studied by UV ab- 
sorption, and which are substantially weaker than the 
forbidden lines, remained under the detection threshold 
of our spectra. This limits our ability to make a direct 
comparison of abundances determined from recombina- 
tion lines for our PNe. However, although our spectra 
did not detect recombination lines originating from any 
of the ions in Tables [7] and [9] from which UV resonance 
absorption was observed, we did observe C II, N II, and 
O II recombination lines. Any information that can be 
obtained from an analysis of these recombination lines is 
potentially useful. Of the above ions only the O II re- 
combination lines originate from a parent ion for which 
forbidden lines were observed, viz., +2 , so the column 
densities inferred from these lines are considered in the 
following section. 

5.3. Recombination Line Column Densities 

The results of Table [9] show that the absorption and 
forbidden emission column densities agree within the un- 
certainties of measurement error, inaccuracies in the val- 
ues of temperature and density, and inhomogcneities that 
cause the adjacent lines of sight to sample different com- 
ponents of the nebulae. This agreement indicates that 
nebular analyses based upon forbidden emission lines 



TABLE 10 
Recombination Line Abundances 



Line(Multiplct) 



(erg cm -2 



ADF 



- 1 ) 



NGC 6543 



O 11 A4069.62,.88 (10) 


5.52 


-14) 


2.24 


A4072.15 (10) 


4.32 


-14) 


1.88 


A4110.79 (20) 


1.16 


-14) 


5.01 


A4153.30 (19) 


2.72 


-14) 


3.61 


A4317.14 (2) 


2.10 


-14) 


2.59 


A4345.56 (2) 


2.47 


-14) 


3.28 


A4349.43 (2) 


2.85 


-14) 


1.57 


A4638.86 (1) 


2.91 


-14) 


2.59 


A4641.81 (1) 


7.34 


-14) 


2.87 


A4649.14 (1) 


8.59 


-14) 


1.75 


A4650.84 (1) 


3.27 


-14) 


3.21 


A4661.63 (1) 


3.20 


-14) 


2.47 


Average 






2.8±0.9 


Tc 1 


O 11 A4069.62,.88 (10) 


4.75 


-16) 


1.97 


A4072.15 (10) 


4.26 


-16) 


1.82 


A4153.30 (19) 


1.58 


-16) 


2.12 


A4317.14 (2) 


2.11 


-16) 


2.88 


A4345.56 (2) 


2.93 


-16) 


4.09 


A4349.43 (2) 


2.98 


-16) 


1.66 


A4641.81 (1) 


7.08 


-16) 


2.72 


A4650.84 (1) 


3.73 


-16) 


3.48 


A4661.63 (1) 


4.37( 


-16) 


3.18 


Average 






2.5±0.9 



a Number in parentheses in an exponent. 

yield heavy element abundances that are the same as 
those derived from absorption lines-the key result from 
this study. That said, what conclusions can be drawn 
from our sample of PNe about nebular abundances based 
on recombination lines? Do our objects show the same 
discrepancies exhibited by other PNe? 

We have detected a number of the same O II recom- 
bination lines from NGC 6543 and Tc 1 that have been 
studied extensively in PNe over the past decade and used 
to determine relative Q +2 abund ances (|Liu et alj 120001 : 
iRobertson-Tessi fe GarnetH [2005) . Since the [O III] for- 
bidden lines are strong in both objects it is straightfor- 
ward to determine the abundance of +2 as derived from 
the two types of lines. No recombination lines were ob- 
served in the spectrum of the very low ionization nebula 
He2-138. 

The permitted O II lines from multiplets 1, 2, 10, 19, 
and 20 have been shown to be populated by recombi- 
nation and to yield, among themselves, consistent val- 
ues of the Q +2 abundance in H II regions and PNe 
(jTsamis et alj|2003t [Wesson. Liu, fe Barlowl 120051) . We 
have used the extinction-corrected intensities of these 
lines, the observed intensities of which are shown in Ta- 
ble Uni to compute the +2 abundance relative to that 
determined from the [O III] AA5007, 4959 lines. This ra- 
tio is referred to as the "abundance discrepancy factor" 
(ADF), and for PNe and H II regions always has values 
that are greater than unity. Using the same cross sections 
and p roced ures described by Roberts on- Tessi fe Garnettl 
(|2005h and I Wesson. Liu, fe Barlowl (|2005f ). we have de- 
termined ADF values from the individual O II lines, 
and these are listed in the last column of Table [lOl 
The resulting mean values of the ADFs for NGC 6543 
and Tc 1 are 2.8 and 2.5 (0.45 and 0.40 dex), re- 
spectively. The mean ADF of 2.8 found here for 
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NGC 6543 is consistent with the previous determinations 
of ADF = 3.0 an d 3.8 from other lines of sight through 
this same nebula (iKingsburgh, Lopez. fc PeimbertllT996l ; 
IWesson fe Liul l2004h . Thus, at least two of our objects 
show the same discrepancies between the recombination 
and forbidden line abundances for +2 that are typical 
of PNe, and we have not by chance studied nebulae for 
which the ADFs are close to unity. 

Since there is good agreement between the absorp- 
tion and forbidden line column densities in our objects, 
one can infer from the above results that recombination 
lines are likely to produce emission column densities that 
are significantly higher than those derived from absorp- 
tion lines. The final column of Table [5] shows that the 
mean of the forbidden emission line column densities is 
marginally larger than that of the absorption column 
densities for each of the PNe. The recombination line col- 
umn densities would produce a larger discrepancy. How- 
ever, because the ADFs for NGC 6543 and Tc 1 are of 
the same size as the uncertainties in the column den- 
sities, a larger sample of PNe is needed, especially in- 
cluding some objects with relatively large ADFs, before 
a definitive statement can be made that recombination 
abundances are not correct. Given that we do not de- 
tect any recombination lines from parent ions for which 
we measured UV absorption lines, a direct comparison of 
absorption and recombination line column densities for 
the same ions is likely to remain elusive. Realistically, the 
only ions in Tables [7] and [9] that are likely to be parent 
ions of detectable recombination lines are S + and S +2 . 
With deeper spectra it should be possible to observe the 
S II recombination lines in our PNe, however the rele- 
vant recombination coefficients are not known and are 
very difficult to calculate with any accuracy (P. Storey, 
private communication) . 

Emission line analysis of a large number of PNe has 
shown that recombination lines originating from C +2 , 
N +2 , +2 , and Ne +2 ions all yield roughly the same 
relative abundances among themselves as do the col- 
lisionally excited forbidden and intercombination lines 
from these same ions, and that the recombination 
lines consistently indi cate higher abundances with re- 
spect to H and He fe obcrtson-Tcssi & Gar nettl 120051 : 
IWesson. Liu, fc Barlowfl2005t iLiu et al.ll2006h . The dis- 
crepancies do not appear to arise from problems with 
the atomic data, rather they seem to be linked to char- 
acteristics that are specific to the nebulae. Thus, the 
ADFs for doubly ionized CNONe tend to be approxi- 
mately the same in individual objects, and they vary in 
step with each other from one nebula to the next al- 
though there are exceptions to this rule, e.g., NGC 6720 
(jGarnett fc Dinersteinll2001h . 

If the agreement among the CNONe ADFs also ap- 
plied to S +2 one could use the ADFs we have derived 
from the O II lines in NGC 6543 and Tc 1 to infer the 
recombination line column density for S +2 , based on the 
[S III] emission measure. However, the ADFs for ele- 
ments in the 3rd row of the periodic table, including sul- 
fur, are virtually unknown because so few recombination 
lines ar e detected and the r elevant cross sections are not 
known (jBarlow et al.ll2003h . The only 3rd row ion for 
which a recombination abundance has been determined 
is Mg +2 , from Mg II l ines h aving been measured in ten 
PNe bv lBarlow etafl (|2003). They found the Mg +2 /H+ 



abundances for their objects to show little evidence for 
enhancement over the solar value, contrary to the large 
+2 /H + enhancements derived from the O II recombi- 
nation lines in the same PNe. They conclude that the re- 
combination line abundance discrepancies may be a phe- 
nomenon restricted to ions of the 2nd row of the periodic 
table, viz., C, N, O, and Ne, and not exhibited by 3rd 
row ions. 

Our current study shows nonetheless that the electron 
densities and temperatures deduced from the usual for- 
bidden line analysis are correct over a range of ionization 
zones that should include the regions where C +2 , N +2 , 
+2 , and Ne +2 reside. This means that the large ADFs 
for the second row elements cannot be reflecting errors in 
the forbidden line abundances due to the use of incorrect 
values of n e or T e . It would be necessary to find another 
mechanism that would affect the forbidden line abun- 
dances derived for C +2 , N +2 , +2 , and Ne +2 , but not 
those found for S + and S +2 . In our opinion this makes 
factors affecting the recombination line abundances al- 
most certainly the cause of the abundance discrepancies. 

6. SUMMARY 

The results reported here are derived from a limited 
sample of PNe and are based upon observations that 
may not be extended in the near future unless STIS is 
repaired and put back into service on HST. For this sam- 
ple we find that the forbidden lines yield absolute abun- 
dances for C I, Fe II, Ni II, O I, S II, and S III that are 
consistent with those derived from their absorption lines 
along adjacent sight lines. Within the uncertainties in 
the line intensities and calculations, the good agreement 
between the column densities derived from the forbidden 
emission lines and the UV absorption lines for the three 
PNe represents a validation of both types of analysis. It 
strengthens confidence in the abundances derived from 
forbidden emission lines in spite of discrepancies with the 
abundances derived from high level permitted recombi- 
nation lines from the same ions, and is the primary result 
of this investigation. 

Although recombination lines were detected in only 
two of the three objects in this study, those two PNe 
do show factor 2.5-2.8 discrepancies between the +2 
abundances derived from forbidden lines and those from 
recombination lines. This demonstrates that we have 
studied PNe in which the abundance discrepancy prob- 
lem exists. Not being able to independently measure an 
abundance for +2 from its UV resonance lines, which 
fall outside of the HST wavelength range, we cannot 
confirm that the recombination abundances for +2 are 
anomalously high. Nor can we use the similarity in the 
abundance discrepancy factors of C +2 , N +2 , and Ne +2 
with that of +2 to make a comparison of the inferred 
recombination abundances of these ions with those from 
an absorption analysis because their UV resonance lines 
also fall outside of the HST wavelength range. The one 
doubly ionized ion for which we do have good absorption 
data, S +2 , did not have recombination lines detected in 
our nebular spectra. 

We have shown that the electron densities and temper- 
atures deduced for a wide range of ionization levels do 
give correct abundances using the forbidden lines from 
ions of other elements within the same ionization zone as 
0+ 2 (and C+ 2 , N+ 2 , and Ne+ 2 ). In particular, the for- 
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bidden line abundances for S +2 are in good agreement 
with the absorption line abundances for S +2 . This is an 
important constraint since any explanation of the ADF 
for +2 and other second row elements that implicates 
errors in the forbidden line abundances would have to 
invoke a mechanism that does not affect the forbidden 
line abundances of other ions such as S +2 in the same 
ionization zone. 
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